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Abstract 

This paper describes a dynamic model that was developed to predict changes in turbine tip 
clearance - the radial distance between the end of a turbine blade and the abradable tip seal. The clearance 
is estimated by using a first principles approach to model the thermal and mechanical effects of engine 
operating conditions on the turbine subcomponents. These effects are summed to determine the resulting 
clearance. The model is demonstrated via a ground idle to maximum power transient and a lapse-rate 
takeoff transient. Results show the model demonstrates the expected “pinch point” behavior. The paper 
concludes by identifying knowledge gaps and suggesting additional research to improve the model. 


Nomenclature 


Symbol 

Units 

Description 

A 

m 2 

area 

E 

Pa 

Young’s modulus 

F 

N 

force 

L 

m 

unstressed blade length 

N 

RPM 

rotational speed 

P 

Pa 

pressure 

T 

°C 

temperature 

V 

m 3 

volume 

c 

J/kg-°C 

specific heat 

h 

W/m 2 -°C 

convection heat transfer coefficient 

k 

W/m-°C 

thermal conductivity 

1 

m 

stressed blade length 

q 

W 

heat transfer rate 

r 

m 

radius 

t 

sec 

time 

u 

m 

deformation 

w 

m 

width 

a 

1/°C 

thermal expansion coefficient 

S 


tip clearance 

P 

kg/m 3 

density 

V 


Poisson’s ratio 

CO 

rad/sec 

angular speed 

<J 

Pa 

stress 

£ 

m/m 

strain 

n 


film cooling effectiveness 

71 


ratio of the perimeter of a circle to its diameter 
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1. Introduction 

Turbine blade tip clearance continues to be a concern in the design and control of gas turbines. 
Ever increasing demands for improved efficiency and higher operating temperatures require more 
stringent tolerances on turbine tip clearance. A successfully implemented active tip clearance control is 
expected to have benefits of increased efficiency, reduced specific fuel consumption, and additional 
service life. 

The focus of this study is to develop a simplified model from first principles that grossly captures 
the dynamic change in turbine tip clearance during take-off. Such a model can aid in identifying and 
understanding the primary mechanisms involved in the turbine tip clearance problem and their inherent 
frequency response. This work serves to identify knowledge gaps that must be addressed in order to 
facilitate design and development of an actuator and control strategy for active clearance control. 
Furthermore, this work is a step towards a higher fidelity model suitable for virtual prototyping * of a 
controller/actuation mechanism. The purpose of beginning with this first-principles approach is to 
eliminate the necessity for a finite- element model that would require detailed information about geometry 
specific to a particular engine system. The model developed herein is intended to be more generic. 

The main objectives of this study can be summarized as follows: 

1. Identify and characterize the primary mechanisms involved in turbine tip clearance dynamics; 

2. Identify the known and u nk nown parameters necessary to predict turbine tip clearance; and 

3. Develop a simplified model that can be used to: 

(a) reasonably represent turbine tip clearance during common gas turbine transients (take- 
off, reburst, etc.), 

(b) determine actuator performance requirements for active clearance control, and 

(c) serve as a baseline for virtual prototyping of a controller and actuation system. 

Previous NASA efforts in collaboration with industry have quantified tip clearance transients 

using both experimental and numerical results [1, 2, 3]. One example is the Energy Efficient Engine 
studies [1] conducted in the late 1970’s and early 1980’s. Tinder these studies, detailed analyses and 
experiments were performed to identify the dynamic characteristics of mechanisms impacting turbine tip 
clearance. Due to proprietary issues, these reports, however, do not describe analysis used to derive 
numerical results. Much of this work has been reserved for classified reports which document analyses 
performed using proprietary finite- element tools, and few publications document the combined effects of 
both thermal and mechanical stresses experienced by the rotor, blade, and shroud. Little has been 
published on the stress analysis used to predict deflections of the rotor, blade, or shroud. Most of the 
finite-element studies in the area focus on long term effects like creep and crack propagation and do not 
focus on prediction of deflections due to thermal and rotational transients during flight. Significant 
experimental and numerical work has been published on the film cooling of turbine blades [4, 5, 6, 7, 8]. 
For NASA to conduct fundamental research that will facilitate future technologies in active turbine tip 
clearance control, a valid model or a test rig is necessary. The test rig can be potentially costly to build 
and run. A model-based simulation is less costly, and is useful in designing the test rig and associated 
experiments. If sufficiently accurate, it can be used to virtually prototype the controller/actuator 
mechanism. 

This document is divided into four chapters, introduction, model development, implementation 
and results, and discussion. The section on model development is presented in three sub-sections based 
on the three system elements composing the turbine - the shroud, rotor, and blade. Each element is 
presented by describing the parameters and boundary conditions necessary for analysis and by detailing 
assumptions used to facilitate results. Analyses used to determine thermal and mechanical deflections are 
described for each element. This is followed by a section in which parameter values are defined for 
simulation purposes; and preliminary results of clearance predictions are presented and discussed. The 


* Virtual prototyping refers to the process by which models and simulations are employed to virtually 
design and test mechanisms and/or controls. 
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paper concludes with a discussion in which the knowledge gaps of most concern are identified and 
summarized, and suggestions for further model development are presented. 

2. Model Development 

As shown in Figure 1, the proposed model incorporates three basic elements - the shroud, turbine 
rotor (or disk), and turbine blade. In order to predict deflections of each element due to thermal and 
mechanical stresses the temperature, pressure, and force distributions in each element must be modeled. 
The block diagram shows the parameters necessary to calculate these temperature and stress distributions. 
An engine model, not discussed in this report, is employed to provide speed, temperature, and pressure 
transients for the shroud, rotor, and blade sub-models. Each sub-model predicts deflections due to 
thermal and mechanical stresses. As shown in Equation 1, the relative change in the time -varying 
geometry of each sub-model is then used to calculate the overall change in tip clearance. 

S(t) = r shroud ( t ) - \r mtor (I) + l hlade (f )] 

= ( r a + u s i + u s2 ) - [(r 0 + u r[ + u r2 ) + {L + u bl + u h2 )\ 

Elere, r sfiroud (t), r rotor (t), and IhimieU) are, respectively, the shroud inner radius, rotor outer radius, and blade 
length as a function of time. Note that, r a , r 0 , and L are the initial geometric state of the shroud, rotor, and 
blade respectively; while the subscripted w’s denote deformations due to thermal and mechanical forces. 


Shroud 

Parameters: 

• geometry of unstressed shroud (w, r a , r b , and r c ) 

• material properties of thermal barrier (k s , c a , and p„) 

• material properties of shroud (at,, k b , c b , and pb) 

• convection coefficient of thermal barrier (b a ) 

• convection coefficient of shroud (h b ) 

Inputs: 

• turbine inlet airflow temperature (Tiurbine) 

• turbine inlet airflow pressure (P turbine ) 

• compressor discharge temperature (7 com p res sor) 

• compressor discharge pressure ( P compressor) 
Outputs: 

• shroud/thermal barrier inner radius (r sbr0 ud) 

• average shroud temperature (T S hroud) 


Turbine Engine Transient 

Outputs: 

• compressor/turbine angular speed (w) 

• compressor discharge pressure (P amr , ressor ) 

• compressor exit temperature ( Tcompressor) 

• turbine inlet pressure (P tU rbine) 

• turbine inlet temperature (r turb i„ e ) 


> 


Turbine Rotor 

Parameters: 

• unstressed rotor geometry (r) 

• material properties (a, c, E, k, p and v) 

• rotor convection coefficient (h) 

Inputs: 

• compressor discharge temperature (T compressor ) 

• angular speed (cj) 

Outputs: 

• rotor tip radius (r rotor ) 

• average rotor temperature (7 r otor) 


r shroud 


r" rotor 


Clearance 

i 

tbl a de 


► 


Turbine Blades 


Parameters: 


• geometry of unstressed blade (L and A) 


• material properties (a, c, E, k, p, and v) 


• convection coefficient for film cooling (h) 


• film cooling effectiveness (p) 


Inputs: 


• temperature of primary airflow (r tur bine) 


• temperature of film cooling airflow (r com prsssor) 


• angular speed ( cj ) 


Outputs: 


• stressed blade length (/ b t a de) 


• blade temperature (r b i ade ) 



Figure 1. Block diagram for tip clearance model 
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Equations for calculating these deflections will be given in the subsections that follow. 

2.1. Shroud 

As shown in Figure 2, it is assumed that the general geometry of the shroud is a ring-like 
structure where the inner surface is coated with an abradable material that also functions as a thermal 
barrier. However, the shroud is not a single continuous ring-like structure. Instead it is made up of a 
series of arcs joined together to form a ring, which may negate the usefulness of this assumption. This 
issue shall be addressed in § 2.1.4 with an alternative formulation of the shroud growth due to thermal 
stresses on the casing to which the shroud is attached. The inner surface of the abradable material at 
radius r a is exposed to heated gas at a temperature approximately equal to the turbine inlet temperature, 

T turb ine • The outer surface of the shroud at radius r h is exposed to compressor discharge air which is 
assumed to be at the compressor discharge temperature, T compressor . As shown in Figure 2, pressure 
distributions on the inner and outer surface of the shroud are adapted from Fattime and Steinetz [9]. For 
simplicity, pressure in the tip clearance region is assumed to vary linearly between the turbine inlet and 
exit pressures. Compressor bleed air is used to purge the space between the shroud and the casing. To 
facilitate results, it is assumed that the temperature difference between the compressor discharge and the 
shroud outer surface is negligible. 

For the purposes of this study, the abradable material is presumed to function primarily as a 
thermal barrier and not as a structural member. It is assumed, therefore, that the abradable material acts 
as a coating and does not impose a significant stress on the inner surface of the superalloy layer as the 
shroud grows due to thermal loads. Additionally, because the shroud is attached to the case in the form of 
a series of arcs, it is assumed that the abradable material layer maintains its approximate thickness even 
when thermally stressed (i.e. the relative distance between the bond radius and inner wall, r c -r a , is 


p 

1 Purge 
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assumed to remain relatively constant). Due to its conductive properties, the abradable layer dissipates 
much of the turbine heat, thus reducing the temperature otherwise experienced by the superalloy surface 
at radius r c . This temperature, T c , is the thermal boundary condition on the inner surface of the alloy used 
to calculate the deflection due to thermal stresses. The heat transfer at the inner surface of the abradable 
material and outer surface of the alloy are assumed to be uniform ( quasi-steady) to facilitate a simple 
model. Heat is convected to the inner surface of the abradable material at wall temperature T a with a heat 
transfer coefficient of h„. The convection coefficient at the outer surface of the shroud is h and the 
associated surface temperature is T h . 


2.1.1. Shroud Heat Transfer 

To facilitate a closed form solution for the heat transfer across the shroud, the responses near the 
shroud surfaces are assumed to be governed by a semi-infinite formulation of the transient heat 
conduction equation where the temperature at the surface, T w , is given by the classical solution [6, 7]: 


T -I 

w i 

T-T 


= 1 


■exp 


Jr 

pck 


-t 


erfc 


h 


sfpck 


t 


( 2 ) 


where T r and 7) are the reference temperature and initial wall temperature, respectively. This formulation 
is appropriate as long as the transient temperature penetration does not exceed the thickness of the wall 
material. For a single-flow system, the reference temperature is the temperature of the mainstream flow 
which for the inner surface of the shroud is assumed to be T turbine and for the outer surface T compressor . This 
formulation is used to calculate the inner and outer shroud wall temperatures, T a and T b , respectively. 

Using a thermal resistance formulation [10] as shown in Figure 2, one can estimate the bond 
temperature, T c , at the interface between the abradable material and the superalloy: 


T a ~T b _ T a -T c 
In (r c /r a ) | In (r b /r c ) In (r c /r a ) 
2nk a w 2 7tk b w l7ik a w 


T-T 

1 a fb 

I , In (r b /r c ) k a ' 

In (r c /r a )k h 


(3) 


Additionally, by representing the superalloy layer of the shroud as a hollow cylinder and 
assuming one-dimensional heat conduction in the radial direction under quasi-steady state conditions, the 
heat diffusion equation simplifies to 


1 d (, dT \ 

kr — 

r dr ^ dr 


which, when integrated and combined with the boundary conditions T(r c )=T c and T(r h )=T h , results in the 
radial temperature distribution 


T{r) = T c +{T h -T c ) 


In (r/r c ) 
In (r b /r c )' 


(4) 


2.1.2. Shroud Deflection due to Thermal Stresses, u s i 


To enable prediction of shroud deflection due to thermal stresses, it is assumed that the stresses 
and displacement do not vary over the width of the shroud. The formulation presented herein is 
equivalent to that given in [1 1] for a thin circular disk with a hollow center (like a ring). According to 
Timoshe nk o [11], the radial and tangential stresses (cr. and <J 0 ) must satisfy the equation of equilibrium 


d& r | Jf ® 8 

dr r 


= 0 


(5) 


*The application of the semi- infinite solution should be justified in future work and is used here primarily 
to facilitate initial results. 
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The strain components in the radial and tangential directions are 

„ 1 . . du 

£ r =ocT + -((J,-vo R ) = — 

E dr 


and 


1 


u 


e g = oO' +—(<7 g -v<7 r ) = - 

E r 

The above strain equations can be solved for the radial and tangential stresses: 

E 

I O J.1/L 

■'d 


and 


1-V 

E 


[e, ,+V£ g ~(\ + v)aT] 
\£ e + V£ r -{\ + v)al] 


Given the above stress equations and noting that £ r = du/dr and £g= u/r the equilibrium equation 
(Equation 5) may be written as 


d_ 

dr 


1 d(ru) 
r dr 


n , dT 
= (1 + v)a—. 
dr 


c. 


( 6 ) 


Integrating twice yields the radial deflection 

1 r' 

u{r) = (1 + V)a— Trdr + C x r H — - 

y J r c y 

where r c is the lower limit of the radius or in this case the unstressed inner wall radius of the shroud 
superalloy layer. Assuming no thermal stress at the inner and outer boundaries of the superalloy portion 
of the shroud (surfaces stresses due to pressure will be considered separately in § 2.1.3) and recalling the 
radial temperature distribution (Equation 4), C x and C 2 are: 


and 


Ci 


c 2 = 


a(\-v) 


T + 


1 


\ r t 


■ K 2ln(/;/r) 


(T b ~T c ) 


c> J 


a(\ + v)r; 


T + 


1 


K r b ~ r c 


2 In (r b /r c ) 


(T b -T c ) 


c> J 


(7) 


( 8 ) 


Because the abradable material is thought to serve mainly as a thermal barrier, the shroud 
deflection is primarily due to the superalloy layer and thus should be measured from the unstressed bond 
radius, r c . The thickness of the abradable material layer can later be subtracted since, as previously 
mentioned, it is assumed to remain relatively constant. Thus the shroud radial deflection due to thermal 
stresses is 


»vi = u(r c ) = a r c 


T + 


1 




■r 2 In (r h /r c ) 


(T h -T c ) 


c> ) 


(9) 


2.1.3. Shroud Deflection due to Pressure Differential, u S 2 

To predict deflection of the shroud at the bond radius, r c , due to pressure differential, 
Timoshenko’s hollow cylinder equation is adapted [11]: 

1 


E(r b 2 ~r c 2 ) 


[d -y)r e (Pf 


P 0 r;)-(l + V)r c r b fP 0 -P i 


)]• 


( 10 ) 


NASA/TM— 2003-212226 


6 



In Equation 10, P, is the mean internal pressure and P„ is the outer surface pressure. Per Figure 2, they 
are equivalent to (P in + P oul )/ 2 and P respectively. 


2.1.4. Alternative Shroud Growth Formulation 


Alternatively, the deflection of the shroud can be thought of as limited by the casing growth. The 
shroud can be then modeled as a series of arcs rigidly attached to the casing as shown in Figure 3. The 
casing grows due to thermal and pressure differentials. Similar to the previous formulation, the surface 
temperatures on the outer and inner surfaces of the casing, T e and T d , can be calculated from mainstream 
flow temperatures and experimentally or numerically determined convection coefficients using 
Equation 2. The outer casing wall, r e , is assumed to be exposed to ambient temperature and pressure (T„ 
and FT). If the casing is grossly modeled as a hoop-like structure then deflection due to thermal 
differential can be derived as was done in § 2.1.2. The deflection due to thermal stresses would, 
alternatively, be as shown in Equation 1 1 . 


’.si 


ar, 


T d + 


1 


r e ~r d 2 In (rjr d ) 


(T e ~T d ) 


d jj 


( 11 ) 


Similarly, the alternative deflection due to pressure differential can be derived as 

1 


*s2 


E(r e 


2 \ 
■ r d ) 


[(\-V)r d (Prj -P o r 2 )-(\ + v)r d r 2 (P -P t \ . 


( 12 ) 


where P, is now the pressure inside the cavity, P purge , and P„ is the ambient pressure, P 0 «. 

In this model, a higher bandwidth mechanical clearance control system is presumed to replace the 
slower thermal clearance control systems found on many current commercial engines. Therefore, active 
cooling of the case is not considered. The deformation of the structure that connects the shroud to the 
case is another effect that is not developed here. The structure deforms as the differential pressure across 
the shroud changes, and so, should be included in the alternative formulation. 


T, & T T„&P„ 
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2.2. Rotor 


Detailed studies of rotor-stator systems have been conducted by Owen [12], However, in the 
attempt to arrive at a simplified model representing the gross clearance dynamics, this initial effort has 
focused on modeling the first-order effects. For this initial study, the rotor structure is modeled as a 
rotating disk of uniform thickness, and only deflections due to changes in centrifugal forces and purge air 
temperature are considered. Compressor discharge air is generally used to purge and cool the cavity 
around the rotor. Only the small tip surface of the rotor is exposed to the high temperature core flow. 
Thermal expansion is assumed to be dominated by heat transfer between the rotor disk and the 
compressor discharge air, at temperature T compressor . Thus, as suggested by Figure 4, the rotor is assumed 
to be insulated to heat transfer from the turbine blades simplifying the analysis. Compressor discharge air 
flows over the majority of the rotor’s surface area. 

2.2.1. Rotor Deflection due to Thermal Stresses, u r i 

As with the outer shroud surface temperature, the rotor surface temperature is assumed to be 
governed by the semi-infinite formulation (i.e., Equation 2) where the reference temperature is the 
compressor discharge temperature. Because most of the rotor’s surface is exposed to the discharge 
temperature, the temperature through the width of the rotor is assumed to be relatively constant across the 
surface and approximately equal to the surface temperature (i.e., T rotor =T compressor ). Note that in this 
simplified model, windage heating, which may produce a significant temperature gradient between the 
hub and the tip, is ignored. 

As with the outer shroud surface, the heat transfer convection coefficient is not known, but it is 
assumed to be similar in magnitude to that given for the inner shroud surface. Since the rotor temperature 
is assumed to be relatively constant throughout, the thermal strain analysis is significantly simplified: 

du u 

£ r = (TT = — and £ e = Off - - . 

dr r 

Therefore, the approximate rotor deflection due to thermal stresses is 

u rl =aTr 0 (13) 

where a r and r 0 are the rotor thermal expansion coefficient and unstressed radius, respectively. 
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2.2.2. Rotor Deflection due to Centrifugal Forces, u ,2 

To estimate the rotor tip radial deflection due to centrifugal forces, a simplified form of the 
rotating flat disk equation from reference 1 1 is used: 

W '-2 =4f( l ~ V r)pyr 0 3 (14) 

4 E 

where V r , p„ and is,- are material properties and (O=(2n/60)N is the rotor angular speed in rad/sec. 

2.3. Blades 

Of the three basic components, the most documented is the turbine blade. The blades are stressed 
by both centrifugal forces and thermal expansion. Because of the high temperatures within the turbine 
stage the turbine blades are film cooled by compressor discharge air. As shown in Figure 5, compressor 
discharge air that flows through the blade exits film cooling holes on the leading edge of the blade and is 
introduced into the boundary layer at the surface of the blade. The blades are of length L and cross- 
sectional area A c . 

2.3.1. Film Cooling of Turbine Blades 

Film cooling is employed to reduce the heat transfer from hot air in the turbine at the turbine inlet 
temperature, T turbine , to prevent surface temperatures in excess of the blade material’s melting point. In 
film cooling, a secondary flow is introduced into the boundary layer at the surface of the blade [13]. The 
temperature of the film cooling fluid, 7}, is much lower than the mainstream temperature, T m . The blade 
surface temperature, T b i ade , is determined using a semi-infinite formulation similar to Equation 2 with the 
exception that the reference temperature, T t , is not the mainstream temperature as in two temperature 
flows. The reference temperature for film cooling is generally unknown and depends on the supply 
temperatures of the two interacting flows and the degree of mixing that occurs [6, 7]. In such cases, the 
reference temperature and convection coefficient, h, must be determined experimentally or numerically. 
Using experimental data, a film cooling effectiveness, similar to that shown in Equation 15, is determined 
in terms of the two flow temperatures, T, n and 7} [7]. 



compressor 


▲ 

L 

T 


Figure 5. Blade schematic 
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(15) 


71 = 


T -T 

r m 

T -T 

1 f 1 m 


In the case of the turbine blades, the mainstream fluid is at T turhine and the film cooling flow is 
assumed to be at the compressor discharge temperature, T compressor . 


2.3.2. Blade Deflection due to Thermal Stresses, Uhi 

The blade is relatively thin, and thus the thermal gradient through the thickness of the blade is 
assumed negligible. Therefore, the blade is assumed to have a uniform temperature equal to the surface 
temperature, T h i ade . This simplifies the thermal stress analysis and results in a simple linear relation: 

u ht = aTL . (16) 


2.3.3. Blade Deflection due to Centrifugal Forces, Ub 2 


Like the rotor, the blade is stressed by centrifugal forces. To simplify results, the maximum blade 
and rotor growth are thought to be sufficiently small to assume that the centrifugal force, F c , is primarily a 
function of the rotor’s angular velocity, oft), and does not vary significantly due to the rotor and blade 
growth. The deflection can hence be determined from a simple stress strain relation: 


£ 


a _ F c t A c _ m{L + r 0 )co 2 _ pL(L + r 0 )co 2 
E~ E ~ EA C ~ E 


U b2 


_ AL _ PL 2 {L + rW 

E 


A L 

L 


(17) 


2.4. Clearance Calculation 

Having obtained expressions for the time varying deflections of the various turbine sub-models 
due to thermal and mechanical forces, these can now be summed to obtain the resulting tip clearance. 
Equation 1, reproduced here for convenience, shows how the relative change in the time-varying 
geometry of each sub-model is used to calculate the overall change in tip clearance. 

S(t) = r shmud ( t ) - \r rotor (0 + l hlade (f )] 

= (fa + U sl + U sl) ~ [(>0 + U , \ + 2 ) + (L + U bl + U b2 )\ 

Here, S(t) is the overall tip clearance as a function of time. Variables r shroud , r wtor , and l blade are, 
respectively, the shroud inner radius, rotor outer radius, and blade length as a function of time. Note that, 
r a , r 0 , and L are the initial geometric state of the shroud, rotor, and blade, respectively; while the 
subscripted w’s denote time-dependent results of previously described deformation calculations. 

3. Model Implementation and Results 

In this section specific parameters used in implementing a dynamic simulation containing the 
previously presented tip clearance model are discussed, along with results from that simulation. Details 
of the actual implementation are not given here, but are given in the Appendix included at the end of this 
report. 


3.1. Implementation 

In this section, values are given for the various parameters required to implement the tip clearance 
model described in this paper. The parameters include geometry, material properties, and heat transfer 
coefficients. A consistent set of data containing all of the necessary model parameters was unavailable. 
Therefore, it was necessary to use estimates for many of the parameters. What follows is a discussion of 
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the sources, confidence level, and associated assumptions used in determining the various parameter 
values. 


3.1.1. Shroud 

Pressure distributions on the inner and outer surface of the shroud as shown in Figure 2 are 
adapted from Lattime and Steinetz [9]. For simplicity, pressure in the tip clearance region is assumed to 
vary linearly between the turbine inlet and exit pressures. Per Lattime and Steinetz [9], these pressures 
are approximately 70% and 30% of compressor discharge, respectively. The resulting average pressure, 
0.5 P compressor, is used on the inner surface of the shroud to determine the shroud deflection due to pressure 
differential. Compressor bleed air is also used to purge the space between the shroud and the casing. The 
resultant pressure on the outer shroud surface is assumed to have an even distribution approximately 80% 
of P compressor- In this preliminary model, delays between the pressure signal at the compressor exit and the 
subsequent signal in the turbine are assumed negligible. 

The shroud parameters and material properties are given in Table 1. Subscript a denotes 
materials and surface conditions of the abradable material, while subscript b denotes the same for the 
superalloy. For the purposes of this model the superalloy layer is assumed to be made of Inconel 718 and 
the abradable material of a partially stabilized Zirconium Oxide. Parameter values are rated on a level of 
confidence from 1 to 5; 5 indicating near certainty. Experimental approximations of the convection 
coefficient on the inner surface of the shroud (i.e., in the tip clearance region) are given in [8] and match 
reported experimental results for the blade tip [4, 5, 6, 7]. The convection heat transfer coefficient for the 
outer surface of the shroud is not known at this time but to facilitate results is assumed to be similar in 
magnitude to that on the inner shroud surface. 

For the alternative shroud growth formulation, the inner casing wall, r d , is exposed to the 
compressor discharge temperature, T compressor , and approximately 80% of the compressor discharge 
pressure, 0.8P compressor . Convection coefficients for the case are not readily available. 


Table 1. Shroud parameters 


Parameter 

Description 

Value 

Confidence 

r a 

thermal barrier inner radius 

38 cm (~15 in) 

4 

r b 

shroud outer radius 

38.6 cm (-15.2 in) 

3 

r c 

shroud inner (bond) radius 

48.3 cm (-15.1 in) 

3 

w 

shroud width 

5.1 cm (-2 in) 

4 

Pa 

abradable material density 

5030 kg/m 3 

3 

k a 

abradable material thermal conductivity 

2.93 W/m-°C 

3 

Ca 

abradable material specific heat 

456 J/Kg-°C 

3 

Pb 

shroud density 

8220 kg/m 3 

5 

kb 

shroud thermal conductivity 

24.9 W/m-°C 

5 

Cb 

shroud specific heat 

645 J/Kg-°C 

5 

a h 

shroud thermal expansion coefficient at 538°C 

1 5. lxl 0" 6 1/°C 

5 

ha 

convection coefficient at inner surface 

500 W/m 2 -°C 

4 

h b 

convection coefficient at outer surface 

1000 W/m 2 -°C 

2 


3.1.2. Rotor 

Rotor disks in most modem gas turbines are fabricated using newer members of the Rene family 
of superalloys. Some rotors [1] are fabricated from Rene 95 which has properties similar to Inconel 718. 
For this work, properties of Inconel 718 are used. Table 2 gives parameters used in the rotor sub-model. 
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Table 2. Rotor parameters 


Parameter 

Description 

Value 

Confidence 

ro 

unstressed radius 

32.9 cm (-13 in) 

3 

w r 

width 

5.1 cm (-2 in) 

3 

Pi 

density 

8220 kg/m 3 

4 

E 

modulus of elasticity 

139 GPa 

4 

K 

Poisson’s ratio 

0.4 

4 

k r 

thermal conductivity 

24.9 W/m-°C 

4 

c r 

specific heat 

645 J/Kg-°C 

4 

a r 

thermal expansion coefficient 

1 5. lxl 0‘ 6 1/°C 

4 

hr 

convection coefficient at surface 

500 W/m 2 -°C 

2 


3.1.3. Blades 

The blade parameters, material properties, and parameter confidence levels are given in Table 3. 
For this initial model, blades are assumed to be made of a material similar to Inconel 718. Turbine blades 
found in modem gas turbine engines are generally fabricated as single crystals using the Rene family of 
superalloys. The incoiporation into the clearance model of properties for these superalloys would be 
appropriate for future studies. 

Experimental values for both the convection coefficient and effectiveness are presented in [6, 7] 
and numerical estimations of the convection coefficient can be found in [14, 15]. For the purpose of this 
study, a bulk or blade-surface-averaged effectiveness and convection coefficient are used. Given the 
experimentally determined convection coefficient and effectiveness. Equations 2 and 15 can be used to 
determine the blade surface temperature, T b i ade . 


Table 3. Blade parameters 


Parameter 

Description 

Value 

Confidence 

L 

unstressed length 

5 cm (-1.97 in) 

3 

A c 

cross-sectional area 

6 cm 2 (-0.93 in 2 ) 

3 

P 

density 

8220 kg/m 3 

5 

E 

modulus of elasticity 

139 GPa 

5 

V 

Poisson’s ratio 

0.4 

4 

k 

thermal conductivity 

24.9 W/m-°C 

5 

c 

specific heat 

645 J/Kg-°C 

5 

a 

thermal expansion coefficient at 
871 °C 

16.4x1 O' 6 1/°C 

5 

h 

convection coefficient at surface 

1000 W/m 2 -°C 

3 

il 

film cooling effectiveness 

0.5 

3 


3.2. Results 

In this section the transient clearance is predicted for a ground idle-to-maximum power transient 
and for a lapse-rate take-off transient. As shown in the upper plot of Figure 6, the ground-idle to 
maximum power transient consists of a step in the throttle setting between the two operating conditions. 
A lapse-rate take-off transient involves ramping the throttle setting from ground-idle to about 90% of 
maximum power, performance of final flight checks, and then ramping up to maximum power. The 
throttle transient for the lapse-rate takeoff is shown in the upper plot of Figure 8. 

For the ground idle-to-maximum power transient, the time history of the engine speed is shown 
on the upper axes of Figure 6. Associated temperatures transients for the high pressure turbine and the 
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Figure 6. Deflections and clearance for transient from ground idle (G/I) to maximum power 
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high pressure compressor are shown on the upper axes of Figure 7. Similar plots for the lapse-rate-take- 
off transient are shown on the upper axes of Figure 8 and Figure 9. Recall from Figure 1, that these 
parameters serve as inputs to the shroud, rotor, and blade sub-models. It is important to note that a 
transient simulation of a large commercial engine was not available. Instead, a simulation of a small 
military engine was used to generate the input transients. The transients were then scaled (as shown) to 
approximate the dynamics of a large commercial engine, for which the clearance model was developed. 


In addition to the transient temperature inputs, Figure 7, shows the shroud temperatures and the 
basic element temperatures. Note that the shroud temperature, T shroud , plotted on the lower axis, is just the 
average of the bond and outer wall temperatures, T c and T b , respectively. As evidenced by the figure, the 
dynamic temperatures generally follow the expected trend - a relatively slow response compared to the 
input temperature transients. The outer wall temperature approaches an expected steady-state value 
between 1200-1300 °F. However, the bond temperature does not appear to approach an expected value of 
near 1800 °F. This is likely due to three primary reasons: 

1. the convection coefficient for the outer shroud surface is not kn own or well documented so the 
assumed value may be incorrect, 

2. the actual properties of the abradable material are not known. A change in the thermal 
conductivity of the abradable material significantly affects the bond temperature, and 

3. the input engine transients may not match those for a larger commercial engine. 

Data on the lower axes also show that the blade wall temperature, T b!ade , stays well below Inconel 71 8’s 
melting point of 2370 °F. 

Figure 6 shows the individual deflections and overall clearance with reference to the engine speed 
transient. As expected the rotor initially responds more quickly due to the centrifugal forces induced by 
the engine speed transient. The shroud catches up and grows more rapidly due to thermal stresses than 
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Figure 8. Deflections and clearance for lapse-rate takeoff transient 
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Figure 9. Temperatures for lapse-rate takeoff transient 


either the rotor or blades. According to this formulation, the blades contribute the least to the clearance 
change. Like the rotor, they grow initially quickly due to centrifugal forces, but their growth is 
significantly less than the rotor which should be expected considering the relatively short length of the 
blades compared to the rotor radius. The predicted clearance generally follows the expected trend for a 
take-off-transient. The clearance is expected to reduce about 30 mils and the recovery due to the thermal 
growth of the shroud is expected to be about 10 mils. The model predicts a reduction of approximately 20 
mils and a recovery of nearly 5 mils. Exact ranges for a given engine would depend on the parameters 
and transient conditions specific to that engine. Nonetheless, the results are promising. Similar 
comments can be made of results garnered for the lapse-rate-take-off transient (refer to Figure 8 and 
Figure 9). 


4. Discussion 

The goal of this effort is to develop a relatively simple physics-based model with sufficient 
accuracy for control development research. The formulation presented here, though grossly simplified, 
appears to reasonably predict the expected change in turbine tip clearance during a take-off-transient. 
Moreover, the predicted temperatures and deflection respond, for the most part, as physically expected. 

Discrepancies do exist, but most deviations from expected response can likely be attributed to 
generalized or inaccurate information for many of the parameters used in this model. Studies by Owen 
[12] provide detailed information on various mechanisms effecting tip clearance, particularly with respect 
to the rotor and shroud. It is not the purpose of this initial model to incorporate all of the details 
contained therein. However, these studies provide a basis for incoiporating additional details that could 
improve the accuracy of clearance calculation. 
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The true shroud thermo-mechanical response is particularly unclear, especially where interactions 
with the case are concerned. The turbine blades are relatively straight-forward, and, in terms of heat 
transfer, have been studied extensively. The blade analysis presented is, therefore, expected to be more 
accurate than that of the shroud or rotor. The rotor can be further complicated due to its rather non-trivial 
geometry and the fact that, in some gas turbines, the rotor is cooled by air at different temperatures on 
either of its sides. 

In light of this study and the questions raised, the following tasks must be addressed to further the 
application and usefulness of this physics-based modeling approach: 

• Confirm materials and associated properties. Because materials are application dependent, the 
material properties of Inconel 718 have been used to approximate those of various superalloys 
used in the manufacture turbine components. A target engine/turbine should be identified and 
accurate material properties used to further develop and validate this model. In particular, the 
abradable material must be more closely researched. The temperature dependence of many of 
the material properties should also be considered in future models. 

• Determine unknown convection coefficients. The outer shroud surface and rotor surfaces are 
cooled by compressor discharge air and the convection process is not well documented. 

• Determine role of interaction between casing and shroud. The shroud growth is likely limited 
by the casing. Documentation on casing heat transfer is not readily available. 

• Quantify roles of windage heating and heat transfer at rotor tip play in determining rotor- 
growth. It was assumed that the rotor temperature is primarily determined by heat transfer from 
cooling air discharged from the compressor. However, the effects of windage heating should be 
investigated as it can cause significant variation in the radial temperature distribution of the 
rotor. Significant heat may also be transferred from the heated blades to the circumferential tip 
surface of the rotor. 

• Identify thermal and mechanical response of abradable material. The abradable material was 
assumed to function mostly as a thermal barrier. Its exact material properties are not known. It 
was assumed that the material did not respond like an alloy resulting in a thermal expansion. 
This assumption must be critically examined. 

Future work should involve more detailed research into those parameters and mechanisms least 
documented. The current model should also be refined and its fidelity quantified in order to assess its 
usefulness for the design of a tip clearance control system. This assessment would be most meaningful if 
real engine data were used as a basis for comparison. 
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Appendix 


This appendix includes the Matlab© m-file and a series of block diagrams from the Simulink© model 
used to create the plots presented in § 3. As shown below, the m-file, GI2MaxPower . m, contains all 
parameters required for the clearance model used to produce the results. Most parameters are 
documented within the comments of the m-file. Input transients used in the calculations were loaded 
from file, GI2MaxPower .mat, into the Matlab© workspace. The MAT-file includes engine speed, 
temperature, and pressure transients used to produce the results shown in Figure 6 and Figure 7. To 
generate the results shown in Figure 8 and, a different MAT-file containing input transients for a lapse- 
rate takeoff was used. The input transients, not given here, were scaled in an attempt to make them more 
representative of a larger gas turbine. The Simulink© models were organized to resemble the block 
diagram shown in Figure 1. 

GI2MaxPower. nt: Matlab © m-file 

% GI2MaxPower . m (created 071102, updated 011603/plotting) 

% Turbine tip clearance parameters 

% Parameter file for calculation of expansion due to 
% thermal and mechanical stresses. 

clear 

help GI2MaxPower 

% Load MAT-files 
%load eng3 

load GI2MaxPower %Ground Idle to Max Power Transient 

%Define Geometry, Thermal & Mechanical Properties, 

% Turbine dynamic inputs 

RtoK = 0.555556; % R to K, temperature conversion 

time = eng.time*2; % slow down transient input 

RPM = eng.xnh*2.45 - 21500; % increase engine speed range 

t41_K = eng.t41*RtoK*0 . 85; % K, turbine inlet temperature 

% needed to scale compressor temperature to make sure it matched 

% published temperatures for cooling air around shroud (1250 F) 

t3_K = eng . t3*RtoK* 1 . 25 ; % K, compressor discharge temperature 

t41_C = t41_K - 273.15; % C, turbine inlet temperature 

t3_C = t3_K - 273.15; % C, turbine inlet temperature 

Po = eng . p3* 68 94 . 7 6*0 . 8 ; % Pa, compressor discharge pressure 

Pi = (0 . 5/0 . 7) *eng.p4*6894 . 76; % Pa, average turbine stage pressure 


P3 = eng . p3* 68 94 . 7 6 ; 

"o 

Pa, compressor exit pressure 


% Unstressed parameters 

rO = 32.9/100; 

"6 

m, rotor tip radius 


L = 0.050; 

"6 

m, blade unstressed length 


Ac = 6e-4; 

"o 

m^2, blade cross-sectional area 

P = 0.110; 

"6 

m, blade perimeter 


R0 = r0+L+0 .0014; 

% Properties of Inconnel 

718 

% m, thermal barrier inner radius 

rho = 8 . 22/1000*100 A 3; 

"o 

kg/m~3, density 


nu = 0.4; 

"o 

poissons ratio 


c = 645; 

'o 

J/Kg-C, specific heat 


alphal = 15.1e-6; 

"o 

1/K, thermal expansion coef. 

at 538 

alpha2 = 16.4e-6; 

"6 

1/K, thermal expansion coef. 

at 871 
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k = 24.9; % W/m-K, thermal conductivity at 871 C 

E = 139e9; % Pa, modulus of elasticity at 870 C 

2'2-S'9-S'9'S'S'2'2-2'9-2'2-S'2'S'2-9- 
ooooooooooooooooooo 

% Shroud parameters 
w = 2*0.0254; 
a = R0; 

b = a + 0.25*0.0254; 
rc = (a+b) /2; 

Aa = 2*pi*a*w; 

Ab = 2*pi*b*w; 

VI = pi*rc A 2*w - pi*a A 2*w; 

V2 = pi*b A 2*w - pi*rc A 2*w; 
ha = 500; 
hb = 1000; 

% Zirconium Oxide 
rhoa = 5.03/1000*100 A 3; 
ca = 0.109*4184; 
ka = 2.93; 
alphaa = 8.2e-6; 

% Inconel 718 
rhob = rho; 
cb = c; 
kb = k; 

alphab = alphal; 

% Initial conditions 
TaO = 1 4 1_C ( 1 ) ; 

TbO = 1 3_C ( 1 ) ; 

ooooooooooooooooooo 

% Blade parameters 
V = Ac*L; 
m = rho*V; 

As = L*P ; 
h = 1000; 
eta = 0.5; 
alpha = alpha2; 

% Initial conditions 
wO = RPM(l) *2*pi/60; 

Ti_blade = (1 - eta)*t41 

2'9-2'9-2'9-S'2-2'9-2'9-2'9-2'9-2'2-2- 
ooooooooooooooooooo 

% rotor parameters 


hr = 500; 

"5 

W/m A 2-K, convection coefficient 

wr = w; 

"6 

m, width 

Ar = 2*pi*rO A 2; 

"o 

m A 2, convection surfe area 

Vr = Ar*wr; 

"6 

m A 3, volume 

rhor = rho; 

"o 

kg/m A 3, density 

Er = E; 

"o 

pa, modulus of elasticity 

nur = nu; 

"o 

poisson's ratio 

alphar = alphal; 

"6 

1/K, thermal expansion coef. 

cr = c; 

"6 

J/Kg-C, specific heat 

kr = k; 

"o 

W/m-K, thermal conductivity 

% Initial Condition 
TrotorO = t3_C(l); 

"o 

C, initial rotor temperature 
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% m A 3, volume 
% kg, mass 

% m A 2, surface temperature 
% W/m A 2-K, film convection coefficient 
% bulk film cooling effectivness 
% 1/K, thermal expansion coef. 

C (1) + eta*t3_C (1) ; 


% shroud with 

% m, thermal barrier inner radius 
% m, shroud outer radius 
% m, shroud inner radius 

% m A 2, thermal barrier inner surface area 
% m A 2, shroud outer surface area 
% m~3, thermal barrier volume 
% m^3, shroud volume 

% W/m A 2-K, conv. coef. thermal barrier 
% W/m A 2-K, conv. coef. shroud 

% kg/m A 3 
% J/kg-K 
% W/m-K 

% 1/K, thermal expansion coef. at 1000 C 

% kg/m A 3, density 
% J/Kg-C, specific heat 
% W/m-K, conductivity at 871 C 
% 1/K, thermal expansion coef. at 538 C 



, 100 ) ; 


% run the simulation 

[T,X,Y] = sim ( ' ClearanceModel9 ' 


RPM 
TempC 
Def IMm 
ClearMm 
ShTempC 


Y(:,l); 

(Y ( : , 2 : 6) -32) *5/ 9; 

Y ( : , 7 : 9) *0 . 0254; 

Y ( : , 10) *0 . 0254; 

(Y ( : , 11 : 15) -32) *5/9; 


%temperature in degrees Celcius 
%deflections in millimeters 
%clearance in millimeters 
%shroud temperatures in degrees Celcius 


S'2'2'2-2'2'2'2'2'2'2'2'2'2'2'2'2'2'2- 

ooooooooooooooooooo 

% plot the results 

set ( 0 , 1 Def aultLineLineWidth ',1.5) 

time = [eng. time; 100]; PLA = [eng. pc; 50]; 

figure ( 1 ) ; elf 

PLApcnt = 100* (PLA-21) /29; 

subplot (3,1,1); [ha, hi , h2 ] =plotyy (T, RPM, time, PLApcnt ) ; 
set (ha, 'Xlim' , [0 100] ) ; %Common x-axis 
axes (ha ( 1 ) ) ; 

set (gca, ' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 

set (gca, ' Ylim' , [7000 17000], ' Ytick' , [7000:2500:17000]); 

ylabel ( ' N (rpm) ' ) ; 

hl= legend ( ' N (rpm) ' , ' Throttle ' , 0 ) ; 
hlc=get (hi, ' children ' ) ; 
axes (hi) ; ylim=get (hi, 'Ylim' ) ; 

line (get (hlc (2 ) , ' Xdata ' ) , get (hlc (2 ) , ' Ydata ')-0.45*diff (ylim) , 
'Color ' , get (h2, 'Color'), 'Line Style', get (h2, 'Line Style' ) ) ; 
axes (ha (2 ) ) ; 

Xlim = get (gca, 'Xlim' ) ; 

Ylim = [0 100]; YTick = [0 100]; YTlbl = [ ' G/I ' ; 'Max ' ] ; 
set (gca, ' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 
set (gca, ' Ylim ', Ylim, ' YTick ' , YTick, ' YTickLabel ', YTlbl ) ; 
ylabel (' Throttle Setting', 'Rotation', - 
90 , ' VerticalAlignment ' , ' bottom ' ) ; 

set (h2 , ' LineStyle ' , '--'); 
subplot (3,1,2); 

%Left Y-Axis is metric units 

plot (T, Def IMm ( : , 1 ) , T, Def IMm ( : , 2 ) , ' — ' , T, Def IMm (:, 3) ,'-.') ; 
set (gca, ' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 

Ylim = get (gca, ' Ylim ') ; YTick = Ylim(l) +diff (Ylim) * [0 : 0 .25: 1] ; 

set (gca, ' YTick ' , YTick) ; 

ylabel ( 'Deflection (mm) ' ) ; 

hl2 = legend ( ' Shroud ' , ' Rotor ' , ' Blade ' ) ; 

%Right Y-axis is English Units 
Xlim2 = get (gca, ' Xlim ') ; 

Ylim2 = get (gca, ' Ylim ')*( 1000/25 . 4 ) ; 
ax2 = axes (' position ', get (gca, ' position ')) ; 
h2=plot (Xlim2 , Ylim2 ,'.'); set (h2 ,' visible ',' of f ') ; 

set (ax2 , ' Xlim ' , Xlim2 , ' Ylim ' , Ylim2 , ... 

' YAxisLocation ' , ' right ' , ' color ' , ' none ' , ... 

' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 
ylabel ( 'Deflection (mils ) ' , ' Rotation ' , -90 , ' VerticalAlignment ' , ' bottom ' 
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subplot (3,1,3); 

%Left Y-Axis is metric units 


plot (T, ClearMm) ; 

set (gca, ' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 

Ylim = get (gca, 'Ylim' ) ; 

YTick = Ylim(l) +diff (Ylim) * [0 : 0 .25: 1] ; 
set (gca, ' YTick ', YTick) ; 
ylabel ( ' Clearance (mm) ' ) ; 
xlabel('Time (seconds)'); 

%Right Y-axis is English Units 
Xlim2 = get (gca, ' Xlim ' ) ; 

Ylim2 = get (gca, ' Ylim ')*( 1000/25 . 4 ) ; 
ax2 = axes (' position ', get (gca, ' position ')) ; 
h2=plot (Xlim2 , Ylim2 ,'.'); set (h2 , 'visible', 'off'); 
set (ax2 , ' Xlim ' , Xlim2 , ' Ylim ' , Ylim2 , ... 

' YAxisLocation ' , ' right ' , ' color ' , ' none ' , ... 

' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 
ylabel ( 'Clearance (mils ) ' , ' Rotation ' , -90 , ' VerticalAlignment 
print -depsc2 -tiff GI2MPdef lections ; 


print -dmeta GI2MPdef lections ; 


' bottom ' ) ; 


f igure (gcf +1 ) ; elf 

set (gcf, ' position ', get (gcf-1 ,' position ').*[ 1 . 05 0.96 1 1]); 
subplot (3,1,1) ; 

%Left Y-Axis is metric units 
plot (T, TempC ( : , 1) , T, TempC ( : , 2) , ' — ' ) ; 
set (gca, ' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 

Xlim = [ 0 100]; 

Ylim = [200 1800] ; 

YTick = Ylim ( 1 ) : 400 : Ylim (2 ) ; 

set (gca, ' Xlim ', Xlim, ' Ylim ' , Ylim, ' YTick ', YTick) ; 

ylabel (' Temperature ( A oC) ' ) ; 

legend ( ' T_{turbine } ' , ' T_{ compressor } ' ) ; 

%Right Y-axis is English Units 
Xlim2 = get (gca, ' Xlim ') ; 

Ylim2 = get (gca, ' Ylim ' ) * ( 9/5 ) +32 ; 
ax2 = axes (' position ', get (gca, ' position ')) ; 
h2=plot (Xlim2 , Ylim2 ,'.'); set (h2 , 'visible', 'off'); 
set (ax2 , ' Xlim ' , Xlim2 , ' Ylim ' , Ylim2 , ... 

' YAxisLocation ' , ' right ' , ' color ' , ' none ' , ... 

' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 
ylabel ( ' Temperature ( A oF ) ' , ' Rotation ' , -90 , ' VerticalAlignment ' , ' bottom ' ) ; 
subplot (3,1,2); 

%Left Y-Axis is metric units 

plot (T, ShTempC ( : , 2 ) , T, ShTempC ( : , 3) , ' — ' , T, ShTempC (:, 4 ),'-.') ; 
set (gca, ' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 

Xlim = [ 0 100]; 

Ylim = [200 1400] ; 

YTick = Ylim ( 1 ) : 400 : Ylim (2 ) ; 

set (gca, ' Xlim ', Xlim, ' Ylim ' , Ylim, ' YTick ', YTick) ; 
ylabel (' Temperature ( A oC) ' ) ; 
legend ( ' T_a ' , ' T_c ' , ' T_b ' ) ; 

%Right Y-axis is English Units 
Xlim2 = get (gca, ' Xlim ') ; 

Ylim2 = get (gca, ' Ylim ' ) * ( 9/5 ) +32 ; 

ax2 = axes (' position ', get (gca, ' position ')) ; 

h2=plot (Xlim2 , Ylim2 ,'.'); set (h2 , 'visible', 'off'); 


NASA/TM— 2003-212226 


22 



set (ax2 , ' Xlim ' , Xlim2 , ' Ylim ' , Ylim2 , ... 

' YAxisLocation ' , ' right ' , ' color ' , ' none ' , ... 

' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 
y label ( ' Temperature ( A oF ) ' , ' Rotation ' , -90 , ' VerticalAlignment ' , ' bottom ' ) ; 
subplot (3,1,3); 

%Left Y-Axis is metric units 
plot (T, TempC ( : , 3) , T, TempC ( : , 4 ) , ' — ' , T, TempC (:, 5), 
set (gca, ' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 

Xlim = [ 0 100]; 

Ylim = [200 1000]; 

YTick = Ylim ( 1 ) : 200 : Ylim (2 ) ; 

set (gca, ' Xlim ', Xlim, ' Ylim ' , Ylim, ' YTick ', YTick) ; 
ylabel (' Temperature ( A oC) ' ) ; 
xlabel('Time (seconds)'); 

legend ( ' T_{ shroud] ' , ' T_{ rotor } ' , ' T_{blade } ' ) ; 

%Right Y-axis is English Units 
Xlim2 = get (gca, 'Xlim' ) ; 

Ylim2 = get (gca, ' Ylim ' ) * ( 9/5 ) +32 ; 
ax2 = axes (' position ', get (gca, ' position ')) ; 
h2=plot (Xlim2 , Ylim2 ,'.'); set (h2 , 'visible', 'off'); 
set (ax2 , ' Xlim ' , Xlim2 , ' Ylim ' , Ylim2 , ... 

' YAxisLocation ' , ' right ' , ' color ' , ' none ' , ... 

' xgrid ' , ' of f ' , ' ygrid ' , ' of f ' , ' box ' , ' of f ' ) ; 
ylabel ( ' Temperature ( A oF ) ' , ' Rotation ' , -90 , ' VerticalAlignment ' , ' bottom ' ) ; 
print -depsc2 -tiff GI2MPtemperatures ; 
print -dmeta GI2MPtemperatures ; 

figure (gcf-1 ) ; 
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Figure A-l. Simulink© block diagram of TipClearance Model 
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Figure A-2. Engine sub-model 



Figure A-3. Shroud sub-model 
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Figure A-4. Inner shroud surface temperature, T a , sub-model 
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Figure A-5. Outer shroud surface temperature, T h , sub-model 
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Figure A-6. Rotor sub-model 
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Figure A-7. Rotor surface temperature, T roton sub-model 



Figure A-8. Blade sub-model 
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Figure A-9. Blade surface temperature, T h / ade , sub-model 
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Figure A-10. Deflection calculations sub-model 
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